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Abstract 

We study the gravitational time delay in ray propagation due to rotating masses in 
the linear approximation of general relativity. Simple expressions are given for the 
gravitomagnetic time delay that occurs when rays of radiation cross a slowly rotating 
shell, equation (14), and propagate in the field of a distant rotating source, equation 
(16). Moreover, we calculate the local gravitational time delay in the Godel universe. 
The observational consequences of these results in the case of weak gravitational 
lensing are discussed. 
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1 Introduction 



Consider the gravitational field of astronomical sources in the linear approx- 
imation of general relativity. The spacetime metric may be expressed in the 
quasi-inertial coordinates x M = (ct, x) as g^ = rj^ + h^, where rj^ is the 
Minkowski metric with signature +2. In the absence of perturbing potentials 
h fll/ (x), rays of electromagnetic radiation propagate along straight lines denned 
by dx l jolt = ck l , where k is the constant unit propagation vector of the sig- 
nal. These zeroth-order null geodesies will be employed throughout this work. 
In the exterior gravitational field of astronomical sources, however, the rays 
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bend due to the attraction of gravity as they follow geodesies of the spacetime 
manifold that are null. It turns out that in order to evaluate the gravitational 
time delay to first order in only the null condition is required, i.e. 

da? da? = . (1) 



We are interested in the consequences of this condition for the propagation of a 
ray on a background global inertial frame. Therefore, it follows from equation 
(1) that c 2 dt 2 — \dx.\ 2 = dx^ dx v , where |<ix| 2 = Sij dx l dx j . To first order 
in the perturbation, the bending of the ray may be neglected in evaluating the 
right-hand side of this result, which may be written as c 2 h a p k a k@ dt 2 . Here 
k a = (l,k) is such that rj^ k v = 0. One may therefore write equation (1) 
in the form 

cdt = (l + hi a p k a k^ \dx\ . (2) 



Let the ray propagate from a point P\ : (eti, Xi) to a point P 2 : (ct 2 , x 2 ) in the 
background inertial frame; then, 

i i 7 

t 2 -h = -|x 2 - Xl | + — k a k? / h a p{x) dl , (3) 

C AC J 

Pi 



where dl = \dx\ denotes the Euclidean length element along the straight line 
that joins Pi to P 2 . It follows from equation (3) that the gravitational time 
delay A is given by 

A G = -| h aP {x)k a k p dl . (4) 
Pi 



Note that in this expression h a p may be replaced by h a p = h a p — ^rjaph, since 
Tj a /3k a k 13 = 0. Here h = ti{h^ v ). Moreover, under a gauge transformation 
of the gravitational potential h^ u — * h^ u + e M) „ + e„ )jU , corresponding to an 
infinitesimal coordinate transformation x^ — > x^ — e M , equation (4) simply 
gives the time delay in the new coordinates. 

In physical circumstances where the linearized Einstein equations are appli- 
cable, the time delay Ac may be combined with the proper time r measured 
by an observer moving with velocity v = c/3, 



dr = cdt \]-g m - 2g oi (3 i - gijfrft , 



(5) 



in order to arrive at physically measurable predictions of the theory. 

The gravitational time delay was studied by Shapiro [1] in terms of radar echo 
delay and subsequently by a number of investigators (see [2]- [8] and references 
cited therein). The general derivation of the linear effect (4) given here is free 
of special assumptions and can be employed in most situations of physical 
interest. 

The gravitomagnetic time delay is discussed in sections 2 and 3. The Godel 
universe is considered in section 4; we study the local gravitational time delay 
in this rotating universe. In section 5 we analyze the gravitomagnetic time de- 
lay in different images of the same source due to gravitational lensing. Finally, 
section 6 contains a brief discussion of our results. 



2 Gravitoelectromagnetism 



Let us now consider the exterior field of slowly moving sources such that the 
trace-reversed potential satisfies □ = — (IQirG/c 4 ) T^, once the gauge 
condition h^ v „ = has been imposed. We are interested in the particular 
retarded solution 



h -1^ 



J T ^ Ct J* ' X0 dV , (6) 



which can be used directly in equation (4) to compute the gravitational time 
delay. The generalization of the Shapiro effect to time-dependent situations 
has been the subject of a recent investigation [6]. 

We confine our study here to the stationary field of a slowly rotating astro- 
nomical body such that h^o = 4§ g /c 2 ,h Qi = — 2(A s )j/c 2 and hij = 0(c~ A ), 
where & 9 {x) is the gravitoelectric potential, A g (x) is the gravitomagnetic vec- 
tor potential (V • A g = 0) and we neglect all terms of order c~ 4 including the 
spatial potentials h^. Far from the source 

GM G Jxr 

®g , A 9 — , 7) 

r c r 6 



where M and J are the total mass and angular momentum of the source. 
We define the gravitoelectric field to be E 3 = — V<I> 9 and the gravitomagnetic 
field to be B 9 = V x A g . In developing this analogy with electrodynamics, one 
encounters extra numerical factors that cannot be avoided. For instance, the 
field equation for the standard gravitomagnetic vector potential has an extra 
factor of —4 in comparison with Maxwell's theory [2]. In this connection, we 



use the following convention that is consistent with the gravitational Larmor 
theorem [9] : The standard formulas of classical electrodynamics are applicable, 
except that for a source of inertial mass M, the gravitoelectric charge is M 
and the gravitomagnetic charge is 2M (in units such that G = 1). On the 
other hand, for a test particle of mass m, the gravitoelectric charge is — m 
and the gravitomagnetic charge is —2m. The source and the test particle have 
opposite charges to ensure that gravity is attractive. Moreover, the ratio of 
the gravitomagnetic charge to the gravitoelectric charge is always 2, since 
linearized gravity is a spin-2 field. 

We find from equation (4) that A G = A GE + A GM , where 

A GE = ^ J §gdl (8) 
Pi 



is the Shapiro time delay and 

Pi 

c 3 „ 
Pi 



A GM = -^- / A <T dx (9) 



is the gravitomagnetic time delay, which is simply proportional to the line 
integral of the gravitomagnetic vector potential. 

Let us now suppose that we can arrange via "mirrors" - - these could be 
transponders on spacecraft — to have the rays travel on a closed trajectory 
as in Figure 1. Then the net gravitomagnetic time delay for a closed loop in 
the positive sense (i.e. counterclockwise) is 

&gm = / A g • dx = -|y B g • dS , (10) 



which is simply proportional to the gravitomagnetic flux threaded by a surface 
S whose boundary is the closed loop under consideration Imagine now a 
point P on the trajectory (cf. Figure 1). If the rays travel along the same path 
but in the opposite direction, the net gravitoelectric time delay measured at 
P will be the same while the gravitomagnetic time delay will change sign, 
A GM = —A GM . Thus the total time difference at P for the rays to go around 
the loop in opposite directions is given by 

t+ - t_ = ~ j A g ■ dx = ~ J B g • dS . (11) 

1 It is interesting to note that the direct echo delay (i.e. when the loop degenerates 
to a line) vanishes in the gravitomagnetic case as the effective area is zero. 



A different derivation of this result based on the propagation of electromag- 
netic radiation in the gravitational field of a rotating mass is contained in ref. 
[10] Q 

Let us now estimate this time difference for the case of GPS signals traveling 
around the Earth. We find that t+ - 1_ « -87rGJ e /(c 4 J R e ) « -10~ 16 sec. By 
contrast, the Shapiro time delay amounts to 2 x 10~ 10 sec for clocks in GPS 
orbit [11], [12], so that this gravitoelectric effect may be measurable in the near 
future. In such an experiment, the time of flight of the signal is monitored and 
eventually measured by a control station on the Earth [11], [12]. The corre- 
sponding gravitomagnetic effect is about two million times smaller and can be 
completely neglected at present. Other very small angular momentum effects 
in laboratory-based optical interferometry experiments have been studied in 
[13], [14]. 



3 Gravitomagnetic time delay 



In this section, we provide explicit formulas for Aqm for the case of radiation 
crossing a slowly rotating thin shell of matter as well as propagating in the 
exterior field of a distant rotating mass. 

In electrodynamics, the magnetic vector potential inside a spherical shell of 
uniform charge density with total charge Q and radius Rq rotating with con- 
stant frequency to is given by A = |B x r, where B = ^Qu>/(cR ) is the 
uniform magnetic field inside the shell. The corresponding electric field van- 
ishes. Therefore, the gravitomagnetic vector potential inside a uniform shell 
of mass M, radius Rq and rotational frequency uj is [15], [16] 



2GM 

A a = u) x r , (12) 

9 3cR 



where the center of the shell is the spatial origin of a background inertial frame. 
Thus inside the shell B g = AGMuj j '(3cR ) and E 9 = 0. The gravitomagnetic 
time delay is therefore 

Pi 

AGM t 

A ™ = -3^7 rxdr ■ (13) 

Pi 



2 A sign error in this reference must be corrected: In ref. [10], the relation <I> + — <!>_ = 
uj(t+ — t-) contains an errant minus sign. Therefore in equations (20), (21) and 
related discussion in section 3 of [10], the sign of t + — t_ has to be reversed. 



Let Pi and P 2 be the two points on the spherical shell indicating the points 
at which the ray enters and leaves the sphere, respectively. Then a simple 
calculation shows that 

AGMRq ^ . . . 

A GM = 7^— w • (ri x r 2 ) , (14) 



where f 1 and f 2 are unit vectors indicating the positions of Pi and P2, respec- 
tively, from the center of the sphere. Note that the result vanishes if u>, ri and 
r 2 are in the same plane; in particular, Agm = if the ray passes through the 
center of the spherical shell. 

Let us next consider the exterior field of a rotating mass with its center of 
mass at the origin of spatial coordinates such that r\ and r 2 are both much 
larger than 2GM/c 2 ; then, neglecting higher-order multipole moments we find 
thatQ 

2G F f J (r x dr) 

- . (15) 



V GM 

c 

Pi 
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A detailed, but straightforward, calculation shows that 



2GJ i 1 1\ J • (f 1 x f 2 ) 

A G m = — — + — — — : — — • (16) 

c 4 \T\ r 2 / 1 + ri ■ r 2 



Thus Agm also vanishes in this case if J, and r 2 are in the same plane. 

Let us note that equation (16) is exact and holds in general for the exterior 
of any spherically symmetric slowly rotating mass with constant J. Moreover, 
let d be the impact parameter of the light ray; then, when T\jd » 1 and 
r 2 /d » 1 it is possible to show that Agm — — 4GJ ■ h/(c 4 d), where n is 
a unit vector normal to the plane formed by r x and r 2 and directed along 
i - ! x r 2 . Extending this calculation to the case in which the source is a thin 
rotating spherical shell of constant radius Rq and the light ray penetrates the 
shell such that the portion of the ray inside the shell subtends an angle tp at 
the center of the shell, we find from equations (14) and (16) that for 7*1 >> Rq 
and r 2 >> Rq, A gm ~ — 4GJ • hf(ip)/(c 4 d), where f((p) = tan ( ZE ^ £ ) + \ sin<£> 
varies from /(0) = 1 to f(n) = as (p : —>■ n. The impact parameter of 
the ray is given by d — i? cos f) so that for cp — we recover the previous 
expression for Agm- 

3 The general expression for the gravitomagnetic time delay in the gravitational 
field of higher-spin multipole moments of a stationary rotating body was derived in 
[4], equation (43). A sign error in formula (43) must be corrected: the overall sign 
of this expression must be reversed. 



In two recent papers [17], [18], the gravitomagnetic time delay in different 
images due to gravitational lensing has been studied and a derivation of the 
time delay due to the spin of an external shell has been presented. These 
results, together with a discussion of the separability of time delays of different 
origin in the images, demonstrate that the gravitomagnetic time delay would in 
principle be measurable in the near future. In particular, the gravitomagnetic 
time delay must be taken into account in the analysis of gravitational lensing 
time delay of some extragalactic sources, since the estimated Aqm may exceed 
the present measurement uncertainty of ~ 0.5 day in the lensing time delay; 
we discuss this issue in section 5. Further measurement aspects of A GM are 
explored in [19], [20]. It is hoped that the measurement of A GM might provide 
information about the existence and distribution of rotating dark matter. 

It would be interesting to estimate how much of the gravitational time delay 
between any two local points Pi and P 2 could be due to the cosmological 
content of the universe. To this end, we consider the development of a Fermi 
coordinate system in the neighborhood of a fundamental geodesic observer 
in the standard FLRW model [21]. At an epoch with Hubble "constant" H, 
the gravitational time delay over any local distance r compared to r/c is of 
the order H 2 r 2 /c 2 , which is negligibly small. A similar estimate holds for a 
rotating universe model as explained in the next section. 



4 Godel universe 

The metric of the stationary and spatially homogeneous Godel universe can 
be expressed as [22] 

ds 2 = -dt 2 - 2^2 U(x) dtdy + dx 2 - U 2 (x)dy 2 + dz 2 , (17) 

where U(x) = exp (a/2 Qx) and we use units such that c = 1 unless specified 
otherwise. It turns out that for this metric the Ricci curvature is given by 
R^ u = 2Q 2 u^u U) where = <5 M is the velocity vector of a particle at rest in 
space and coincides with the timelike Killing vector d t . The physical content of 
this model may be thought of as a perfect fluid with velocity and constant 
density and pressure given by p = p = Q 2 /(8nG), where Qd z is the vorticity 
vector associated with the geodesic worldlines of the fluid. Alternatively, the 
universe could be filled with dust of constant density fl 2 /(4:7iG) together with 
a cosmological constant A = — fi 2 [23]. 

To study the influence of this rotating cosmos on the local physics, we establish 
a Fermi normal coordinate system in the neighborhood of a standard observer 
in this model. Imagine therefore an observer at rest in space and comoving 



with the fluid; then, the observer follows a geodesic and its velocity vector 
= A^q-j is the temporal axis of an orthonormal tetrad frame A'V) that is 
parallel transported along the geodesic. In (t, x, y, z) coordinates, such a tetrad 
is given by 



A" (0) = (1,0,0,0) , (18) 

A" (1) = A M (1) cos Vlt + A M (2) sin Vtt , (19) 

A" (2) = -A M (1) sin Vtt + A M (2) cos Vtt , (20) 

A" (3) = (0,0,0,1) , (21) 

where A^ and A^ 2 ^ are defined by 

A% = (0,1,0,0) , (22) 

\»=(-V2,0,U- l (x),0) . (23) 



Here the rotation of the orthonormal triad -V^, i — 1, 2, 3, representing ideal 
gyroscope directions characterizing the spatial Fermi frame, about the z-axis 
with frequency Q has been made explicit. The Fermi coordinate system as- 
signs coordinates X M = (T, X, Y, Z) to an event P in the neighborhood of 
the geodesic worldline under consideration as follows : there exists a unique 
spacelike geodesic connecting P to the worldline at Po such that the two paths 
are orthogonal at Po, i.e. A^ ^ r/^ = 0, where r/^ is the unit tangent vector to 
the spacelike geodesic at Po- Let r be the proper time of the observer at Po, cr 
be the proper length of the spacelike geodesic segment PoP and C, = ^A^ 
be the direction cosines of this segment with respect to the observer's triad at 
Po. Then T = r, X = aC\, Y = 0-C2 and Z = aCz- In these Fermi coordinates, 
the spacetime in the neighborhood of the observer is Minkowskian except for 
the cosmic tidal influence expressed by the gravitational potential F h iiv , i.e. 
9p,v = Vt+v + F h^ v , where 



F h 00 = - F Ro l03 X*Xi , (24) 
2 

hoi — — - RojikX J X , (25) 

F hij = — - F RikjiX k X l . (26) 

Here the Riemann curvature components are the projections of the Riemann 
tensor on the tetrad of the observer 

F Raf3-y8 = R^upa^( a ) A" (/3) A P (7) A CT (<5) . (27) 



It turns out that all of the nonzero components of the Riemann tensor in 
x^ = (t, x, y, z) coordinates can be obtained from 



#0101 = ^ 2 , #0202 = n 2 U 2 (x) , (28) 

#0112 = -V2 Q 2 U(x) , R 12 i2 = 3Q 2 U 2 (x) , (29) 

using the symmetries of the Riemann tensor. It follows that in the Fermi 
coordinates all of the nonzero components are given by 

F #0101 — ^#0202 = ^#1212 = ^ 2 (30) 

via the symmetries of the Riemann tensor. Thus 



F h 00 = -n 2 (X 2 + Y 2 ) , F h u = -\tt 2 Y 2 , (31) 
% 2 = % 1 = W , F h 22 = -hl 2 X 2 , (32) 

are the only nonzero components of the gravitational potential and the metric 
at this linear order is given by 



ds 2 = - [1 + Q 2 (X 2 + Y 2 )\ dT 2 + dX 2 + dY 2 + dZ 2 

-hl 2 (XdY -YdX) 2 . (33) 



It is now straightforward to compute the gravitational time delay (4) using 
the relations (31) - (32). The result is 



a --^{ '"|r" 2) ' 2 + I r -- r ^'" xR -> 2 

+ (QxR 2 ) 2 + (QxR 1 )-(f]xR 2 )]l , (34) 



where ft = QZ and Hi and R 2 are the position vectors of Pi and P 2 , respec- 
tively. 

The relative time delay for a signal from the observer (R x = 0) to a local 
distance R compared to R/c is thus given by (ft x R) 2 /c 2 , which is negligibly 
small on the basis of the current upper limit of Q < 10~ 24 sec -1 on the possible 
rate of rotation of the universe [2], [24]- [2 7]. 



5 Gravitomagnetic time delay and gravitational lensing 

In gravitational lensing, depending upon the relative position and distance of 
the source, observer and deflecting mass, several images of the same source 
may be observed. Numerous examples of gravitational lensing have been dis- 
covered; a well-known case is the gravitational lens Q2237 + 0305, or Einstein 
Cross, where the path of light from a quasar estimated to be at a distance 
of approximately 8 billion light years is bent by the gravitational field of a 
galaxy estimated to be at a distance of about 400 million light years. This 
light bending produces four images of the same quasar as observed from the 
Earth [28], [29]. 

The gravitomagnetic time delay due to the spin of an astrophysical object 
might then be detected in different images of the same source by gravitational 
lensing. Time delay in the images may be generated by the spins of the de- 
flecting object and of other bodies around the path of the light rays, e.g. by 
the spin of an external rotating mass. 

Let us estimate the spin time delay corresponding to some astrophysical con- 
figurations. For simplicity, to derive the order of magnitude of the gravito- 
magnetic time delay we assume that the source, lens and observer are aligned. 
Nevertheless, there is an additional time delay, called geometrical time delay 
[30] , due to the different geometrical paths followed by different rays. Depend- 
ing on the geometry of the system, this additional term may be very large 
and may be the main source of time delay. However, if we compare the time 
delay of photons that follow the same geometrical path we can neglect the 
geometrical time delay, as in the case of two light rays with the same impact 
parameter but on different sides of the deflecting object. For a small deflec- 
tion angle, the contribution to the travel time delay from the deflected path 
length traveled is of the second order in the deflecting potentials; however, 
depending on the geometrical configuration considered, this delay may need 
to be included in the total time delay. Here, for simplicity, we neglect any 
geometrical time delay. Furthermore, there is a relative time delay due to the 
quadrupole moment of the central mass distribution. It is shown in [17] how, 
in the special case of multiple images of the source with the same impact pa- 
rameter, propagating along the same axis, one may in principle have enough 
observables to solve for the angular momentum J and mass M of the central 
deflecting body by eliminating the time delay due to the unknown quadrupole 
moment. Of course, for other configurations in which the source is not exactly 
aligned with the lens and the observer, the difference in the paths traveled 
and the corresponding difference in the Shapiro time delay can be the main 
source of relative time delay; one would then need to model and remove these 
delays between the different images on the basis of the observed geometry of 
the system. 



Let us now calculate the time delay due to the spin of some astrophysical 
sources. For the Sun, GM Q /c 2 ~ 1.477 km and R Q ~ 6.96 x 10 5 km, by 
considering two light rays traveling on the equatorial plane with impact pa- 
rameters b ~ R Q and —b, the relative gravitomagnetic time delay from formula 
(16) is given by A$ M = 8GJ/(c 4 b) ~ 1.54 x KT 11 sec. The time delay due to 
spin of the Sun could then, in principle, be measured using an interferometer 
at a distance of about 8 x 10 10 km by detecting the gravitationally deflected 
photons emitted by a laser on the side of the Sun opposite the detector and 
traveling on opposite sides of the Sun to the interferometer. To derive the time 
delay due to the lensing galaxy of the Einstein Cross [28], [29], we assume a 
simple model for the rotation and shape of the central object. Details about 
this model can be found in [31]. The angular separation between the four 
images is about 0.9", corresponding to a radius of closest approach of about 
650/if 5 1 pc, and the mass inside a shell with this radius is ~ 1.4 x 10 10 h^i M Q 
[29]. Let us assume that GJ/c 3 ~ 10 23 /if 5 2 km 2 , we then have from formula 
(16) that A£|m — 8GJ/(c 4 b) ~ 4 min. Thus, at least in principle, one could 
measure the time delay due to the spin of the lensing galaxy; of course, as in 
the case of the Sun, one should be able to model with sufficient accuracy and 
remove all the other delays due to other physical effects from the observed 
time delays between the images. As a third example we consider the relative 
time delay of photons due to the spin of a typical cluster of galaxies of mass 
M c — 1O 14 M , radius Rq — 5 Mpc and angular velocity cue — 1CT 18 sec -1 ; 
depending on the geometry of the system and on the path followed by the pho- 
tons, we then find relative time delays ranging from a few minutes to several 
days [17]. 

Especially interesting is the case of gravitomagnetic time delay due to the 
rotation of an external sphere. Let us now calculate the time delay in the travel 
time of photons propagating inside a rotating shell corresponding to some 
astrophysical configurations. In the case of the Einstein Cross [28], in order 
to get an order of magnitude of the effect, we assume that the lensing galaxy 
has a radius R ~ 5 kpc; after some calculations based on the model given 
in [31] and integrating formula (14), the relative time delay of two photons 
traveling at distances of 650 pc on opposite sides from the center is given by 
A r Q M ~ 20 min. If the lensing galaxy is inside a rotating cluster, or super- 
cluster of galaxies, to get an order of magnitude of the time delay due to the 
spin of the mass rotating around the deflecting galaxy, we use typical super- 
cluster parameters: total mass M = 10 15 M , radius R = 70 Mpc and angular 
velocity u = 2 x 10~ 18 sec -1 [32]. If the galaxy is in the center of the cluster and 
light rays have impact parameters ~ 15 kpc (of the order of the Milky Way 
radius), the time delay, in the case of constant density p can be determined 
by integrating formula (14) and the result is A t q M ~ 1 day. 

Finally, if the lensing galaxy is not at the center of the cluster of radius R but 
at a distance aR from the center, with < a < 1, by integrating formula (14) 



between the impact parameter r ~ aR and R, we find that when the impact 
parameters r\ and r 2 of two light rays are approximately equal 



Thus, if the lensing galaxy is at a distance of 10 Mpc from the center of the 
cluster, the relative time delay due to the spin of the external rotating mass 
between two photons with T\ — r 2 ~ 30 kpc is Aq 1 m ~ 0.9 day. 

Promising candidates to observe the time delay due to spin are systems of 
the type of the gravitational lens B0218+357 [33], where the separation angle 
between the images is so small, 335 milliarcsec for B0218+357, that the time 
delay between the images is also very small, about 10.5 days for the two images 
of this system. In such cases, the total time delay is comparable, for a favorable 
configuration, to the time delay due to spin. In addition, an Einstein ring is 
observed in the system B02 18+357 and the diameter of the ring is the same 
as the separation of the images. In such a configuration, the Einstein ring can 
provide strong constraints on the mass distribution in the lens; in turn, this 
can be extremely useful in order to distinguish time delays of different origin. 
Since the present measurement uncertainty in the lensing time delay is of the 
order of 0.5 day, the gravitomagnetic time delay might already be observable 
[34]. For example in the system B0218+357 the measured delay is 10.5±0.4 
day. This measurement is possible because the source is a strongly variable 
radio object, thus one can determine the time delay of the variations in the 
images. In the case of B0218+357 it is possible to observe distinct variations 
in the total flux density, percentage polarization and the polarization position 
angle at two frequencies. 



6 Discussion 

Time delay plays an important role in classical wave scattering as well as 
quantum scattering theory (see [35] for a recent comprehensive survey of this 
subject). Here we have confined our attention to the propagation of radiation 
in the JWKB limit, i.e. null rays, in a gravitational field that is considered 
only to linear order. We have then derived expressions for the gravitomagnetic 
time delay of null rays propagating inside a spinning shell and in the field of 
rotating masses. 

Since we have shown that there may be an appreciable time delay due to the 
spin of a body, or external shell, we draw the following conclusions. 

The gravitomagnetic time delay should be taken into account in the modeling 




(35) 



of relative time delays of the images of a source observed at a far point by 
gravitational lensing. This effect is due to the propagation of the photons in 
opposite directions with respect to the sense of the spin of the body, or external 
rotating shell. 

If other time delays can be modeled with sufficient accuracy and removed 
from the observational data, the larger relative delay due to the quadrupole 
moment of the lensing body can be removed, at least in principle and for some 
configurations of the images, by using special combinations of the observables; 
thus, one could directly measure the spin time delay due to the gravitomag- 
netic field of the lensing body. In order to estimate the relevance of the spin 
time delay in some real astrophysical configurations, we have considered some 
possible astrophysical cases. We have studied the relative time delay in the 
gravitational lensing images caused by a typical rotating galaxy or a cluster of 
galaxies. We have also analyzed the relative gravitomagnetic time delay when 
the path of photons is inside a galaxy, a cluster or a super-cluster of galaxies 
rotating around the deflecting body; this effect should be large enough to be 
detectable from the Earth. 

The measurement of the gravitomagnetic time delay due to the angular mo- 
mentum of an external massive rotating sphere might be a further observable 
for the determination of the total mass-energy of the external body, i.e. of 
the dark matter content of galaxies, clusters and super- clusters of galaxies. 
Indeed, by measuring the gravitomagnetic time delay one can determine the 
total angular momentum of the rotating body and thus, by estimating the 
contribution of the visible part, one can determine its dark-matter content. 

These estimates are preliminary because we need to apply the gravitomagnetic 
time delay to some particular, known, gravitational-lensing images. Further- 
more, we need to estimate the size and the possibility of modeling of the other 
sources of time delay. Nevertheless, depending on the geometry of the as- 
trophysical system considered, we conclude that the relative gravitomagnetic 
time delay may be an already observable effect. 
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Fig. 1. Light propagates from satellite 1, is transponded from satellites 2 and 3, and 
returns back to satellite 1 by making a closed counterclockwise loop. Gravitoelectric 
time delay Age and gravitomagnetic time delay Aqm are due to the gravitational 
field of the rotating mass M possessing spin J. Point P is an arbitrary point on the 
loop. The"clock effect" at the point P is given by equation (11) and also discussed 
in [10] using a different theoretical approach. 



